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Abstract 

Let P and S be two disjoint sets of n and m points in the plane, respectively. We 
consider the problem of computing a Steiner tree whose Steiner vertices belong to S, 
in which each point of P is a leaf, and whose longest edge length is minimum. We 
present an algorithm that computes such a tree in 0{{n + m)\ogm) time, improving 
the previously best result by a logarithmic factor. We also prove a matching lower 
bound in the algebraic computation tree model. 

1 Introduction 

Let P and S be two disjoint sets of n and m points in the plane, respectively. A full Steiner 
tree of P with respect to S is a tree T with vertex set P U S' , for some non-empty subset S' 
of S, in which each point of P is a leaf. Such a tree T consists of a skeleton tree, which is 
the part of T that spans S', and external edges, which are the edges of T that are incident 
on the points of P. 

The bottleneck length of a full Steiner tree is defined to be the Euclidean length of a 
longest edge. An optimal bottleneck full Steiner tree is a full Steiner tree whose bottleneck 
length is minimum. In [1] , Abu-Affash shows that such an optimal tree can be computed in 
0((n + m) log 2 m) time. In this paper, we improve the running time by a logarithmic factor 
and prove a matching lower bound. That is, we prove the following result: 

Theorem 1 Let P and S be disjoint sets of n and m points in the plane, respectively. An 
optimal bottleneck full Steiner tree of P with respect to S can be computed in 0((n+m) logm) 
time, which is optimal in the algebraic computation tree model. 
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2 The algorithm 



2.1 Preprocessing 

We compute a Euclidean minimum spanning tree MST(S) of the point set S, which can 
be done in 0(m log m) time. Then we compute the bipartite graph T 6 (P, S) with vertex 
set P U S that is defined as follows: Consider a collection of six cones, each of angle n/3 
and having its apex at the origin, that cover the plane. For each point p of P, translate 
these cones such that their apices are at p. For each of these translated cones C for which 
C fl S 7^ 0, the graph T 6 (P, S) contains one edge connecting p to a nearest neighbor in C fl S. 
(This is a variant of the well-known Yao-graph as introduced in [5].) Using an algorithm of 
Chang et al. [3], together with a point-location data structure, the graph T 6 (P, S) can be 
constructed in 0((n + m) logm) time. 

The entire preprocessing algorithm takes 0((n + m) logm) time. 

2.2 A decision algorithm 

Let A* denote the optimal bottleneck length, i.e., the bottleneck length of an optimal bottle- 
neck full Steiner tree of P with respect to S. 

In this section, we present an algorithm that decides, for any given real number A > 0, 
whether A* < A or A* > A. This algorithm starts by removing from MST(S) all edges having 
length at least A, resulting in a collection Ti,T%, ... of trees. The algorithm then computes 
the set J of all indices j for which the following holds: Each point p of P is connected by 
an edge of T 6 (P, S) to some point s, such that (i) s is a vertex of Tj and (ii) the Euclidean 
distance \ps\ is less than A. As we will prove later, this set J has the property that it is 
non-empty if and only if A* < A. The formal algorithm is given in Figure [TJ 

Observe that, at any moment during the algorithm, the set J has size at most six. 
Therefore, the running time of this algorithm is 0(n + m). 

Before we prove the correctness of the algorithm, we introduce the following notation. 
Let j be an arbitrary element in the output set J of algorithm CompareToOptimal(A). 
It follows from the algorithm that, for each i with 1 < i < n, there exists a point Sj in S 
such that 

• Si is a vertex of 7), 

• (pi, is an edge in T 6 (P, S), and 

• \piSi\ < A. 

We define 7J to be the full Steiner tree with skeleton tree Tj and external edges Si), 
1 < i < n. Observe that, since each edge of Tj has length less than A, the bottleneck length 
of Tj is less than A. Therefore, we have proved the following lemma. 

Lemma 1 Assume that the output J of algorithm CompareToOptimal(A) is non-empty. 
Then X* < A. 



2 



Algorithm CompareToOptimal(A); 

remove from MST(S) all edges having length at least A; 

denote the resulting trees by 7\, T 2 , . . .; 

number the points of P arbitrarily as Pi,P2, ■ ■ ■ ,Pn] 

J:=0; 

for each edge (pi, s) in T 6 (P, S) 

do j : = index such that s is a vertex of 7}; 

if \pis\ < A 

then J := J U {j} 

endif 
endfor; 
for % := 2 to n 
do for each j E J 

do keep(j) := /aZse 

endfor; 

for each edge (p,, s) in T 6 (P, 5) 

do j := index such that s is a vertex of Tj; 

if j G J and < A 

then keep(j) := true 

endif 
endfor; 

J := {j G J : keep(j) = true} 
endfor; 

return the set J 

Figure 1: This algorithm takes as input a real number A and returns a set J. This set J is 
non-empty if and only if A* < A. 
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The following lemma states that the converse of Lemma [T] also holds. 

Lemma 2 Assume that X* < A. Then the output J of algorithm CompareToOptimal(A) 
has the following two properties: 

1. J 7^ and 

2. J contains an element j such that Tj is a full Steiner tree, whose skeleton tree Tj has 
bottleneck length less than X, and in which each external edge has length at most A*. 

Proof. Consider an optimal bottleneck full Steiner tree, let T* be its skeleton tree, and 
denote its external edges by (pi, Sj), 1 < i < n; thus, each Si is a vertex of T*. Each edge of 
this optimal tree has length at most A*. 

We may assume that T* is a subtree of MST(S); see Lemma 2.1 in Abu-Affash [Tj. Since 
each edge of T* has length at most A*, which is less than A, there exists an index j, such that 
T* is a subtree of Tj. We will prove that, at the end of algorithm CompareToOptimal(A), 
j is an element of the set J. 

Let i be any index with 1 < i < n. Recall that the graph Tq(P, S) uses cones of angle 
7r/3. Consider the cone with apex pi that contains Si. This cone contains a point s- of S 
such that (pijS'A is an edge in T e (P,S). (It may happen that s[ = Since \pis[\ < \piSi\, 
we have |sjS^| < \piSi\ < A* < A. 

Consider the path in MST(S) between Sj and s' v It follows from basic properties of 
minimum spanning trees that each edge on this path has length at most \sis\l < A. Therefore, 
s' { is a vertex of the tree Tj. 

It follows from algorithm CompareToOptimal(A) that, when pi is considered, the 
index j is added to J if i — 1, and j stays in J if i > 2. Thus, at the end of the algorithm, 
j is an element of the set J, proving the first claim in the lemma. 

The full Steiner tree Tj, having skeleton tree Tj and external edges (pi, s'^) for 1 < i < n, 
satisfies the second claim in the lemma. I 



2.3 Binary search and completing the algorithm 

Let k denote the number of distinct lengths of the edges of MST(S), and let Ai < A2 < . . . < 
Afc denote the sorted sequence of these edge lengths. Define Ao := and A^+i := 00. Using 
algorithm CompareToOptimal to perform a binary search in the sequence Ao, Ai, . . . , A&+1, 
we obtain an index £ with 1 < £ < k + 1, such that A^_i < A* < Xe. 

Since algorithm CompareToOptimal takes 0(n+m) time, the total time for the binary 
search is 0((n + m) logm). 

Run algorithm CompareToOptimal(A^). Since A* < Xe, it follows from Lemma[2]that, 
at the end of this algorithm, the set J contains an index j such that Tj is a full Steiner tree, 
whose skeleton tree Tj has bottleneck length less than A, and in which each external edge 
has length at most A*. Since Tj is a subtree of MST(S), it follows that each edge of Tj has 
length at most A^_i, which is at most A*. Thus, TJ is a full Steiner tree with bottleneck 
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length at most A*. By definition of A*, it then follows that the bottleneck length of Tj is 
equal to A*. 

Thus, to complete the algorithm, we run algorithm CompareToOptimal(A^) and con- 
sider its output J. For each of the at most six elements j of J, we construct the full Steiner 
tree Tj and compute its bottleneck length A*. For any index j that minimizes A*, Tj is 
an optimal bottleneck full Steiner tree. This final step completes the algorithm and takes 
0(n + m) time. This proves the first part of Theorem [TJ 

3 The lower bound 

In this section, we prove that our algorithm is optimal in the algebraic computation tree 
model; refer to Ben-Or [2] for the definition of this model. 

3.1 The case when n is small as compared to m 

We start by assuming that n = 0(m). We will prove that the problem of computing 
an optimal bottleneck full Steiner tree has a lower bound of f2(mlogm), which is Q((n + 
m) logm). 

Consider a sequence Si, s%, . . . , s m of real numbers. The maximum gap of these numbers 
is the largest distance between any two elements that are consecutive in the sorted order of 
this sequence. Lee and Wu jl] have shown that, in the algebraic computation tree model, 
computing the maximum gap takes fi(mlogm) time. 

Consider the following algorithm that takes as input a sequence Si,s 2 ,...,s m of real 
numbers: 

1. Compute the minimum and maximum elements in the input sequence, compute the 
absolute value A of their difference, and compute the value g = A/(m + 1). 

2. Compute the set S = {(sj, 0) : 1 < i < m}, a set Pi consisting of n/2 points that are 
to the left of (s^, 0) and have distance at most g/2 to (si, 0), a point set P 2 consisting 
of n/2 points that are to the right of (s m , 0) and have distance at most g/2 to (s m , 0). 
Let P be the union of P\ and P 2 - 

3. Compute an optimal bottleneck full Steiner tree T oiP with respect to S, and compute 
the length A* of a longest edge in T . 

4. Return A*. 

Let G be the maximum gap of the sequence si, S2, . . . , s m , and observe that G > g. It is not 
difficult to see that G = A*. Thus, the above algorithm solves the maximum gap problem 
and, therefore, takes Q(mlogm) time. Since n = 0(m), the running time of this algorithm 
is 0(m + n) = 0(m) plus the time needed to compute T ■ It follows that the problem of 
computing an optimal bottleneck full Steiner tree has a lower bound of fi(mlogm). 
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3.2 The case when n is large as compared to m 

We now assume that n = Q(m). We will prove that the problem of computing an optimal 
bottleneck full Steiner tree has a lower bound of fi(nlogm), which is Q((n + m) logm). 

A sequence Pi,P2, ■ ■ ■ ,p n of points in the plane is specified by 2n real numbers. We identify 
such a sequence with the point (p±,P2, ■ ■ ■ ,p n ) in K 2n . For each integer i with 1 < % < m, let 
Q be the point (z, 1). Define the subset V of M 2 ™ as 

V = {(Pl,P2, ■•■,Pn)6R 2 ": {P1,P2, C {Ci,C 2 , . • • ,C m }}. 

For any function / : {1, 2, . . . , n} — )■ {1, 2, . . . , m}, define the point Pf = (c/(i), c/p), . . . , c/( n )). 
Since there are m n such functions /, we obtain m" different points Pf, each one belonging 
to the set V. The set V is in fact equal to the set of these m n points Pf and, therefore, 
V has exactly m" connected components. Thus, by Ben-Or's theorem [2], any algorithm 
that decides whether a given point (pi,p%, ■ ■ ■ ,p n ) belongs to V has worst-case running time 
f2(nlogm). 

Now consider the following algorithm that takes as input a sequence Pi,p%, ■ ■ ■ ,p n of 
points in the plane: 

1. Compute the set S = {(i, 0) : 1 < i < m}. 

2. Let p = (0, 0) and q — (m + 1, 0), and compute the set P' = {p, g} U {pi,P2, ■ ■ ■ ,Pn}- 

3. Compute an optimal bottleneck full Steiner tree T of P' with respect to S. 

4. Set output = true. 

5. For each j with 1 < j < n, do the following: 

(a) Let % be the index such that pj and (z, 0) are connected by an external edge in T. 

(b) If pj 7^ Cj, set output = false. 

6. Return output. 

If the output of the algorithm is true, then each pj is equal to some q and, therefore, the 
point (pi,P2, ■ ■ ■ ,p n ) belongs to the set V. 

Assume that (pi,p 2 , . . . ,p n ) £ V. The (unique) optimal bottleneck full Steiner tree of 
P' with respect to S is the union of (i) the path connecting the points of S sorted from 
left to right (this is the skeleton tree), (ii) the edge connecting p with (1,0) and the edge 
connecting q with (m, 0) (these are external edges), and (iii) edges that connect each point 
Pj of P to the point Cj having the same ^-coordinate as pj (these are also external edges). It 
then follows from the algorithm that the output is true. 

Thus, the algorithm correctly decides whether any given point (pi, P2, ■ ■ ■ ,p n ) belongs to 
V. By the result above, the worst-case running time of this algorithm is Q(nlogm). Since 
m = 0(n), the running time of this algorithm is 0(m + n) = 0(n) plus the time needed to 
compute T . It follows that the problem of computing an optimal bottleneck full Steiner tree 
has a lower bound of Q(nlogm). 

This completes the proof of the lower bound in Theorem [TJ 
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